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A METHOD- OF- MOMENTS  SOLUTION  FOR  DISPERSION  CHARACTERISTICS 
OF  ARBITRARILY-CONFIGURED  TRANSMISSION  MEDIA 

I.  INTRODUCTION 

The  success  of  microstrip  in  microwave  integrated  circuit  appli- 
cations has  caused  considerable  interest  in  the  calculation  of  dispersion 
characteristics  of  these  lines.  A number  of  different  techniques  [l-4] 
have  been  employed  to  obtain  dispersive  effects  of  open  and  shielded 
microstrip-like  transmission  lines  with  rectangular  cross-sections. 

Since  microstrip  becomes  lossy  and  difficult  to  fabricate  at  higher 
microwave  frequencies,  attention  has  focused  on  configuring  new  trans- 
mission media.  In  this  report  we  present  a technique  for  calculating  the 

dispersion  characteristics  of  electromagnetic  wave  propagation  along 

' 

guiding  structures  consisting  of  a finite  number  of  uniform  dielectric 
regions  of  arbitrary  cross-sections  within  a conducting  enclosure. 
Conducting  strips  may  also  be  present  at  the  interface  between  two  di- 
electric regions.  It  is  assumed  that  the  thickness  of  the  conductors 
is  negligible. 

In  Section  II  the  problem  is  formulated  on  the  basis  of  the 
equivalence  principle.  A set  of  linear  integro-dif ferential  operator 
equations  for  the  equivalent  current  sources  are  obtained  by  applying  the 
appropriate  boundary  conditions  at  each  interface.  In  Section  III  the 
operator  equations  are  converted  into  a matrix  formulation  by  the  method 
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of  momenCs  [6,7,8].  The  numerical  methods  used  to  determine  the 
propagation  characteristics  and  results  for  specific  examples  are 
described  in  Section  IV. 

II.  INTEGRO- DIFFERENTIAL  EQUATIONS 

Fig.l  shows  the  generic  cross-section  of  the  guiding  structure 

under  consideration.  N is  the  number  of  discrete  homogeneous, 

isotropic  dielectric  regions  inside  the  conducting  enclosure.  The 

heavy  lines  on  the  interface  between  two  dielectric  regions  denote 

conducting  strips.  The  electric  (E)  and  magnetic  (H)  fields  in  each 

region  will  be  obtained  by  applying  the  principle  of  equivalence  Qs]. 

In  accordance  with  this  principle  the  dielectric  medium  of  the  p-th 

region  (characterized  by  permittivity  « ) is  fictitiously  extended  to 

^ — — 

fill  all  space  and  combinations  of  electric  (J^)  and  magnetic  (M^) 
surface  current  sources  are  conceptually  placed  on  the  boundary  S of 

- P 

the  p-th  region.  Jr  and  M are  to  be  determined  in  such  a way  that 
— * -♦ 

E and  H are  zero  everywhere  outside  and  are  identical  to  the  fields 
Ep  and  Hp  at  each  point  in  the  interior  of  the  p-th  region  for  the 
original  problem  shown  in  Fig.  1,  This  procedure  is  repeated  for  each 
of  the  regions  inside  the  conducting  enclosure.  The  current  sources 
for  the  various  regions  are  not  all  independent  because  of  the  boundary 
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conditions  to  be  satisfied  at  all  the  interfaces.  Fig. 2 symbolically 

shows  the  contour  of  the  p-th  region  and  the  surface  current  distri- 
”*P 

butions  on  it.  M is  taken  to  be  zero  on  the  conducting  segments  of 


the  boundary.  Also  shown  in  Fig. 2 is  a left-handed  coordinate  system 
A A A A 

with  unit  vectors  t , n and  z.  r is  tangential  to  the  contour  c 
P P P P 

(counterclockwise) , n^  inward  drawn  normal  to  the  region  p and  ^ 

perpendicular  to  the  plane  of  Fig. 2.  For  two  adjacent  regions  p and  p', 

we  have  r ■ - t , and  n ■ - ^ . on  the  portion  of  boundary  comnon 
P P P P 

to  both  the  contours  c and  c . . z is  the  same  for  all  regions . 

P P 

^ Q ^ 0 

Since  J and  Mr  are  surface  currents,  we  have  the  relations: 


From  the  equivalence  principle: 


—D  — A — n - 

Jp  (p,  z)  - n X Hp  (p,  z) 


-»tj  -•  A -» n-  -* 

Mp  (p,  z)-  -n  X E p (p,  z) 
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V. 


where  E P and  HP  denote,  respectively,  the  electric  and  magnetic 

fields  at  a point  just  inside  the  contour  cp.  p is  a position  vector 

in  the  tangential  plane.  In  equation  (2)  p is  on  c . 

P 

From  Maxwell's  equations  the  electric  (E1*)  and  magnetic  (H11)  fields 
at  any  point  in  the  crossseccion  of  Fig. 2 may  be  written  as 


-«p  E P«  7X  F P + j epu)AP  +«pV  CDeP  , (3) 

-u  HP«-VXAP+jucoFP+n  V O P , (4) 

o o o in 


where  P ■ 1,  2,  ...,  N . Ap  and  5pgp  are  the  vector  and  scalar  potentials, 
respectively,  due  to  the  electric  current  sources  (J p)  residing  on  the 
contour  cp.  Similarly,  Fp  and  cpmP  are  the  vector  and  scalar  potentials, 
respectively,  due  to  the  magnetic  current  sources  (M  **)  residing  along 
cp.  M>q  is  the  free  space  permeability.  In  equations  (3) -(4)  a time 
dependence  of  e^‘cut  has  been  assumed  for  the  sources. 

The  potentials  A,  F,  CD  and  cp  are  given  by 

e in 


+icut  ( 


Ak  (p,  z,  t)  - HQe 


JP  (p'f  z)C(kP  R)dr'  , 


+iiut 


F (P,  z,  t)  - • e 
P 
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^ M P (p  1 , 


z)G(kP  R)dT' 


(5) 


(6) 


_ icut  ( - _t 

(p,  z,  t)  - -(l/jcu«p)  e $ 7' • Jp  (p  , z , t)G(kPR)dr'  , 


(7) 
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(8) 


(P,  2,  t) 


itut  L _ _ 

- (1/ j(UiQ)e  J>  V'-  MP(p',  z,  t)G(kPR)dT' 


where  R = |p  - p'|  is  the  distance  between  a field  point  ( p ) and  a 

source  point  (p'),  kp  = |u  « cua  + k*|’  and  G(kp  R)  is  the  two- 

pop  z p 

dimensional  Green's  function.  If  (p,oepuj3+  ka)  is  positive,  G(x)  can  be 

(a) 

expressed  in  terms  of  Hq  (x),  the  zero  order  Hankel  function  of  the 


second  kind,  as 


G(x)  - (l/4j)H<a)(x); 


(9) 


while  for  (p.«  ou3  + ka)  negative 
P z 


G(x)  » (1/2TT)  K (x)  , (10) 

J 

where  K^(x)  Is  the  zero  order  modified  Bessel  function. 

The  gradient  operators  v and  v' operate  on  the  field  points  and  the 
source  points  respectively.  The  operator  v'may  be  written  as 


”Vt  + 2 Sz 


. A 3 

+ Z “ - 


(ID 


where  v ^ is  the  tangential  part  of  the  vector.  The  electric  (Ep)  and 

magnetic  (H  ) fields  in  the  p-th  can  be  calculated  by  substituting  Eqs. 

Eqs.  (5)--(10)  in  Eqs.  (3)  and  (4).  We  assume  that  the  sources  vary 
-k  2 

as  e along  the  z-axis.  For  lossless  propagation  k^  = i|3.  By  using 
the  one-dimensional  divergence  theorem  terms  like  (V^  . Jp  )G  in  Eqs. (71 
and  (8)  may  be  expressed  in  the  form  JP  • Vj.  G.  Furthermore,  we  have 
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v ' G(kPR>  = - V G(kPR>  = -R  kP  ’ 
to  to  0 ox 


(12) 


P N - - 

where  x * k^R  and  R * (p  - p')/R  is  the  unit  vector  from  the  source 
point  (p1)  to  the  field  point  (p).  Using  these  relations  and  after 
some  vector  manipulations  the  tangential  components  of  Ep  and  Ip  at  a 
point  P_  just  inside  the  contour  cp  may  be  expressed  in  the  following 
form: 


z o 


■f*1  * 

•if 

[/"■  \ 


< k2  j k kp 
G(x)  h + — /+  -£— ! - 


e k 
p o 


P o 


j IT'  Jj,  r’  . 


R G'  (x) 


A.  A 


dr'  n’-  RG'(x)  MP, 


A A Q 

T . R G (x)  JP 


P 3 

dT'Ji»]^/*  ^ ( G<x)  + r-2-  gxx))  + ^ 

\\  l a / « 1,3 


IdT'^ 


e kaR 
P ° 


A.  A 

• R)(T  • R) 


e k 
P o 


[G"  (x)  - 


kP  R 


G’(x) 


kp  r 

I dr'  G'  (x)  (n 
a J 


• R)””z'  + E*  /< 


dT 1 (n 1 • T)G(x)  Mjj.( 


(13) 


(14) 
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n2  = G*0«0)V> 


'V  ' 


ko  ^ dT'  J^'  n'  * R G,(X’  + ep  I1  + - ,Z  a 1 - dT' 


e k3  I j 

p o / J 


k kP  -i 

ZD  . -v  p \ ■' 

dT 1 M , T 1 • R G 1 (x  I , 


k2  J 


J 


«?- 


, r kP  r -k  f 

= »(  if  I dT'  JZ-  (“  • R)  G'(x)  - if  dT'  J?'  T G(x) 

O J o J 


k k^  r 

dr'^P,  ? - R G*(x) 

ko  J 


+ s ■ dT'  MP  1 (t  • t')[g(x)  + — - — G'(x)J 

r I ' f — . O 


J 


s kaR 
P o 


(k^)2 


e k2 
P o 


'"o'  A A yl  \ r i 

+ 1 — ■“  (T1  • R)  (t  . R)Lg"(x) G”x)] 

. y n ' 


kp  R 


'J 


Here  we  have  replaced  the  variable  M by  M , where 

MP  = -j(e  /ut  f m p 
oo 


so  that  MP  is  of  the  same  dimension  as  J P . 
dielectric  constant  of  the  p-th  region,  k 
wave  vector.  (n,  z)  and  (n1 , t’  3' 


e = s /e  is  the  relative 
p P o 


= cu  e is  the  free  space 
OO 


G(x) 


) denote  unit  vectors  at  the  field 


' 


point  (P)  and  the  source  points  (P'),  respectively.  Also  G' (x)  = 5G(x)/Bx 

and  G"(x)  = BsG/Sx“  with  x = k3  p 
v 0 

If  the  dielectric  medium  of  tue  p-th  region  is  fictitiously  ex- 
tended to  all  space,  then  the  tangential  components  of  the  fields  at  a 
point  P+  just  outside  the  contour  c^  are  given  by  (See  Eq.(2)) 


EP+  = EP”  + j (l  /e  )'2 

Z Z JV00  T 


mP 


e:-  = Er  - j(n  /«  )*  MP 


(18) 


= EP“ 
r r 


o o 


where  E and  H are  given  by  equations  (13)-(19).  The  region  p'  is 
adjacent  to  the  p-th  region  in  the  structure  shown  in  Fig.l.  On  the 
portion  of  the  contour  common  to  regions  p and  p',  /t'  = -9  and  z'  = z. 

In  equations  (20)- (22),  p * 1,  2,  . . . . , N.  Also  for  each  p in  equations 
(22),  p'  = 1,  2,  Np  where  is  the  number  of  regions  adjoining 

the  p-th  region. 

On  substitution  of  equations  (13 ) — (1 9 > in  equations  (20)-(22),  we 
obtain  the  set  of  integral  equations  to  be  satisfied  by  the  current 
sources  JP  , MP  in  all  regions.  JP  and  Mp  in  the  different  regions 
are  not  all  independent  because  of  the  boundary  conditions  in  Equation  (22). 
From  equations  (2)  and  (22)  it  is  seen  that 


-! 


The  propagation  characteristics  are  determined  by  solving  the 


system  of  equations  (13) -(22).  These  equations  are  reduced  to  matrix 
form  by  the  method  of  moments  [6],  The  simplest  approximation  con- 
sists of  using  pulse  functions  for  a basis  and  point  matching  for 


testing.  For  this  purpose  the  contour  c^  is  divided  into  nP  segments 
and  pulse  functions  defined  as: 


fl(p)  * 1 on  g-th  segment, 

P 

5 0 on  all  other  segments. 


The  segments  are  not  necessarily  of  equal  length. 
Let: 


where 


<23) 


(24) 


ny  * number  of  conductor  segments  on  c , 
c p 

n P * nP  - nP  ■ number  of  dielectric  segments  on  c 
a C p 
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The  segment  numbers  have  been  rearranged  such  that  3 = 1,  2,  .....  n^ 


represent  conductor  segments  and  3 3 + 1,  n^  are  dielectric 


segments. 


On  substitution  of  Equation  (24)  in  Equations  (13) -(22)  and 


satisfying  the  resultant  equations  at  the  midpoint  of  each  segment 


we  obtain  the  following  system  of  equations:  Equations  (20)  and  (21) 


become: 


£ 4*5*  £ 


where 


a 3 1,  2,  .....  n1^  , 


p * 1.  2.  .....  N . 


and  equation  (22)  transforms  into: 


R5  6 S + TfiP»  P+nj?  S P+»£ 


* £ < V * £ 


where  the  index  <5  runs  over  all  dielectric  segments  consnon  to  the 


adjoining  regions  p and  p'.  The  index  5'  for  the  contour  cp , refers 


to  the  same  segment  as  the  index  6 for  the  contour  c^.  A set  of 


J 
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equations  similar  to  Equation  (26)  holds  for  each  distinct  pair  (p,  p1) 

of  adjacent  regions.  In  equations  (25)  and  (26),  SP  is  a 1x2  vector 

P 

— P 

representing  the  two  components  of  J , the  electric  current  source  on 

P 

the  0-th  conductor  segment  of  contour  c . The  1x4  vector  f P denotes 

P - P 

-♦  p ^ 

the  two  components  of  J and  two  components  of  , the  current 

sources  on  the  0-th  dielectric  segment  on  contour  c Each  element  Pp 

P Op 

is  a 2x2  submatrix  and  each  element  a 2x4  submatrix.  The  two  rows 

of  these  submatrices  correspond  to  the  two  components  of  Ep  . Similarly 

each  element  Rpg  (and  ) is  a 4x2  submatrix  and  each  element 

T^g  or  TPg  a 4x4  submatrix.  The  four  rows  of  these  submatrices 

-*  n *y 

correspond  to  the  two  components  of  E and  two  components  of  MP  in 

Equation  (22) . The  explicit  forms  of  the  matrix  element  of  Pp  , 

ap 

Qae  ’ Rap*  Ta  0 * *ap  “d  *a0  are  9hown  itl  the  Appendix- 

Equation  (25)  is  a set  of  2nP  equations  in  2nP  + 4np  variables 

c d 

for  the  p-th  region.  Equation  (26)  is  a set  of  4m  equations  (m  being 

the  number  of  dielectric  segments  common  to  regions  p and  p')  in 

JaP  + 4np  + 2np  + 4nP  variables  of  the  regions  p and  p'.  These  equations 

for  all  the  N regions  and  for  all  distinct  pairs  of  adjacent  regions 

3 

can  be  manipulated  in  such  a way  that  the  current  sources  S on  the 
dielectric  segments  are  expressed  in  terms  of  S , the  current  sources 
on  the  conductor  segments  and  then  a matrix  equation  involving  S 
only  is  obtained: 


H S = 0, 
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pvn 


■ - 


mmqmgmga  - 


i 


where  S is  a one-column  vector  with  M components  given  by: 


(28) 


and  H is  a MXM  matrix. 

If  the  cross-section  of  the  guiding  structure  is  symmetrical  about 
an  axis  in  the  transverse  plane,  then  the  problem  can  be  reduced  to 
modes  having  either  even  or  odd  symmetry  with  respect  to  this  axis. 

In  such  a case,  we  have  (for  the  ordering  of  segments  as  in  Fig. 3): 


+ (-l)J+1  SP 

nP  + 1 - 0,  j 


IP  ' •-±(-l)J+1  S 

np  + Ct,j  np  + 1 - a,  j 


(29) 


3P 

nP  +CH,  j+2 


+(-l)J 


3P 

np  + 1 - a,  j+2. 


By  virtue  of  the  definitions  used  for  z and  r',  for  even  symmetry  modes 

-vi  /V 

J and  M remain  unaltered  while  J and  M changes  sign  under  reflection 
Z T T Z 

about  the  symmetry  axis.  The  opposite  is  true  for  odd  symmetry  modes. 

Using  Equation  (29)  the  set  of  equations  given  by  Equations  (25) 
and  (26)  can  be  reduced  to  the  following  form: 


c/2 


P-1 


3 


V'*  > «£.  „P  -p 


-i/a,  n*.  S' 

S-l  r?  +S 

c 


(30) 


RP-  S P + 
~S3  p 


3 -1 


3 -1 


T^e  S ” 

^6,nP+e  Qp  + p 


,P’ 


3-  1 


d n 

RP*  SP’+  > "TP’  D,  SP‘ 

^6'3  3 ^ w 6 ' » np nP'+| 

P-l 


(31) 


where  oc  » 1,  2,  nP/2  and  the  index  & runs  over  the  dielectric 

segments  common  to  regions  p and  p'  ly*  .g  to  the  right  of  the  symmetry 

axis  in  Fig. 3.  The  indices  p,  p'  and  6'  have  the  same  meaning  as  in 

Equations  (25)  and  (26).  In  equations  (30)  and  (31),  the  quantities 

Pp  , RP  and  Rp  are  given  by  the  relation: 

**  CL p •wCt  p 


(-i)J+1  Ap 

(“ 


(32) 
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and  Qr 


have  the  form: 


and  T 


r 

ae 


AP 

~va  , n1'  + f>  I 

K ' ° J» 


1,2 


for  j - 2, 
ij 


( 


The  (+)  and  (-)  signs  in  Equation  (33)  apply  to  the  even  and  odd  symmetry 
modes  respectively.  After  elimination  of  the  variables  S p , an  equation 
similar  to  Equation  (27)  is  obtained  where  H is  now  of  dimension  (M/2)  x 
-(M/2).  Because  of  the  smaller  dimension,  the  computation  time  for  the 
calculation  of  the  determinant  of  the  matrix  H is  greatly  reduced. 


IV.  NUMERICAL  RESULTS 

Since  the  components  of  the  vector  S in  Equation  (27)  are  all 
independent,  a solution  exists  if  and  only  if 


Det  |H(f,  kr)|  - 0, 


where  f » <u/2tt  is  the  operating  frequency  and  k is  the  phase  constant 

in  the  direction  of  propagation.  For  a given  f,  the  propagation  constant 

is  determined  by  finding  k such  that  Equation  (34)  is  satisfied.  There 

z 

will  be  several  values  of  k for  each  f corresponding  to  different  order 

z 

modes.  The  cut-off  frequencies  for  the  different  modes  may  be  obtained 


I 


exact  current  sources.  Due  to  this  approximation,  the  values  of  k^ 
needed  to  satisfy  Equations  (35)  and  (36)  are,  in  general,  slightly 
different.  Therefore,  an  adequate  approximate  solution  is  obtained 
by  requiring  that: 


Det  (H(f , k ))  = minimum.  (37) 

z 

In  actual  numerical  calculation  a few  spurious  roots  occur  in  the  solution 

of  Equation  (34).  The  actual  roots  are  identified  by  the  following  three 

criteria:  (1)  k°  that  gives  the  deepest  local  minimum  in  |Det  h| , 

z 

(2)  real  and  imaginary  parts  of  Det  H should  change  sign  near  k°  , and 

z 

(3)  the  difference  in  the  values  of  k satisfying  Equations  (35)  and  (36) 
is  smallest.  Also,  the  spurious  roots  shift  appreciably  when  the  number 
of  segments  in  a contour  is  changed  Clj. 

To  illustrate  the  accuracy  of  the  present  method,  results  computed 
for  a shielded-microstrip  cross-section  shown  in  Fig. 4 are  plotted  in 
Fig. 5.  The  calculated  dispersion  characteristics  of  the  fundamental  and 
higher  order  even-symmetry  modes  agree  reasonably  well  with  the  theoretical 
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results  in  references  2 and  4.  Although  only  even- symmetry  modes 
are  shown  in  Fig. 5,  the  method  is  applicable  to  all  types  of  modes, 
symmetric  or  otherwise. 


We  next  show  the  dispersion  characteristics  of  a channelized 
suspended  microatrip  [9].  The  cross-section  of  the  structure  is 
shown  in  Fig. 6.  The  channel  located  above  the  conducting  strip  helps 
suppress  higher-order  mode  propagation.  The  structure  has  two 
additional  useful  features:  (1)  reduced  dissipation  loss  LlOJ,  and 
(2)  easier  fabrication  due  to  wider  strip  widths  for  50Q  impedance 
level.  The  calculated  dispersion  curves  are  shown  in  Fig. 7.  The 
lower  three  curves  represent  the  fundamental  mode  for  three 
different  values  of  W/H  (=  1,  3,  4),  the  ratio  of  the  width  of  the 
conducting  strip  and  the  height  of  the  channel  above  it.  The  upper 
two  curves  are  two  higher  order  (even  symmetry)  modes  for  W/H  * 3. 

The  associated  TEM  phase  constants  in  air  and  dielectric  material 
(e^  * 10.0)  are  shown  in  the  figure  for  reference.  As  can  be  seen  in 
Fig.?,  the  phase  constants  for  the  fundamental  mode  at  lower  frequencies 
are  nearer  to  that  in  air  and  at  high  frequencies  go  towards  the  values 
for  the  dielectric  medium.  The  cut-off  frequencies  for  the  next  two 
even-symmetry  higher  order  modes  are  17.2  GHz  and  27.4  GHz,  respectively. 


V.  DISCUSSION 


The  computer-aided  analyses  described  in  this  paper  can  determine 
the  dispersion  and  higher  order  mode  characteristics  for  a wide  variety 
of  transmission  structures  having  different  geometries  and  material 
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parameters.  The  analyses  includes  the  effects  on  propagation  due 
to  an  arbitrarily  shaped  conducting  enclosure.  The  analyses  presented 
here  can  provide  design  information  for  planar  transmission  media 
which  employ  composite  conductor  and/or  dielectric  materials.  The 
analysis  can  be  readily  extended  to  determine  other  propagation 
characteristics  such  as  electric  and  magnetic  field  distributions, 

, impedance  parameters  and  dissipation  losses. 


modal  currents 


APPENDIX 


The  matrices  P,  Q,  R,  R,  T and  T in  Equations  (25)  and  (26)  may 
be  obtained  in  a straightforward  manner  by  substituting  Equations 
(13)- (19 and  Equation  (24)  in  Equations  (20)-(22)  as  stated  in  the 
text. 

The  matrix  elements  P^  for  a + 3 are  given  by: 


where  Rap  - |pQ  - Pg|  and  RQp  - (PQ  - Pp)/*^  • Pa  ^ the  position 

vector  of  the  midpoint  of  the  m-th  segment  and  and  nQ  denote, 
respectively  the  unit  vectors  at  that  point  along  the  tangent  and  the 

normal  to  the  contour.  W_  is  the  width  of  the  m-th  segment.  For  a * S 
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where  Y is  the  Euler  Constant  and  B = -U-if  (e  ka  + k3)  ^0  while  B * 0 
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The  matrix  elements  (P  a are  given  by: 
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for  a j*  0,  and 
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for  OL  - P. 

P 

The  expressions  for  Rag  are: 


( tL  ^ _ f , i - 1,  2 and  j - 3, 
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Fig.  2.  Conceptual  treatment  of  currents  and  fields  for  the  p-th 
region,  r , ft  and  ft  are  unit  vectors  along  three  orthogonal 
coordinates  axes  used  in  the  text,  ft  is  normal  to  the  plane  of  the 
paper  pointing  up. 


Fig.  4.  Cross-section  for  microstrip 
in  a conducting  box 
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